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\q • Abstract 

We address the problem of studying the decay of the infiaton field <f> to another 
scalar field x through parametric resonance in the case of a coupling that involves 
several decay modes. This amounts to the presence of extra harmonic terms in the 
JL | perturbation of the x field dynamics. For the case of two frequencies we compute 

the geometry of the resonance regions, which is significantly altered due to the 
presence of non-cuspidal resonance regions associated to higher harmonics and to 
the emergence of instability 'pockets'. We discuss the effect of this change in the 
efficiency of the energy transfer process for the simplest case of a coupling given by 

l/-, a combination of the two interaction terms of homogeneous degree usually consid- 

ered in the literature. We find that the presence of higher harmonics has limited 
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cosmological implications. 



The success of the inflationary paradigm depends to a great extent on the corre- 
sponding success of the reheating stage that takes place after inflation [H[2] through 
which all elementary particles that exist in the universe were created. During in- 
flation the universe expands exponentially, its matter content is diluted and the 
temperature decreases as the inverse of the exponential of the number of e-folds 
exp N e . Unless some different mechanism is considered for the infiaton decay (as in 
warm inflation [3j, [I]) there has to be some process to raise the temperature to the 
levels required for the nucleosynthesis of the light elements to take place according 
to the standard thermal history of the big-bang universe [S] . 

In most post-inflationary scenarios reheating occurs due to particle production 
by an oscillating scalar field <p. In the simplest models this field is the infiaton field 
that drives inflation. After inflation the scalar field <j> oscillates near a minimum 
of its potential and this triggers a sequence of processes that produces elementary 
particles and eventually restores the temperature J6j \7\ EJ EJ [TQl ttH 12] • 



Since the beginning of the 90's a considerable effort has been devoted to model 
the reheating process |13|, [T^ |6j (a general account can be found in |2j). In most 
models, the first stage of this complicated sequence involves the excitation by para- 
metric resonance of a second scalar field, here denoted x, giving rise to an expo- 
nential increase in the number of x boson particles [15(, [To] (for a comprehensive 
review see also [17]). 

Despite the many contributions regarding the preheating/reheating mechanism 
itself pUl ED E21 E3 EH E51 ESI E7J EE1 [29] and its observational implications [SHED 
[32j[33], some general questions remain to be completely answered. In particular, 
how does the resonant energy transfer depend on the coupling between </> and x an d 
on the inflaton asymptotic dynamics? The latter issue has been studied in depth 
for a large variety of polynomial potentials [El Ell E3 GB] . The former issue is less 
well studied and has been the subject of some recent work JHTJEHl 155] . 

In the simplest of the pre-heating scenarios the inflaton couples to the x field 
through interaction terms of the form h(f>x 2 or g 2 (f) 2 x 2 i that correspond to two 
different decay modes of the scalar field <p into another boson [TJ [SJ. These two 
coupling terms give rise to the same qualitative effects, and are considered as al- 
ternative models. Indeed, in both cases the equation for the scalar field x can be 
reduced to a Mathieu equation, and thus the parametric resonance follows similar 
patterns (even though the numerical values of the model's outcome may be slightly 
different). There are two regimes, a broad resonance regime, in which the ampli- 
tude of the periodic perturbation of the x _ field frequency is of the same order as 
or larger than the frequency of the 4> scalar field, and a narrow resonance regime 
where the amplitude of the perturbation is small. The broad resonance region of 
parameter space includes, for sufficiently large values of the perturbation, the ta- 
chionic resonance regime which has been shown in [37] to be extremely effective in 
transfering most of the energy of the inflation to the x field. 

The main feature that emerges from these studies of pre-heating is that the 
broad and tachyonic resonance regimes gives the predominant contribution to the 
X field energy density. However, there is always the possibility of a contribution 
to the total particle production in the narrow resonance regime when the coupling 
parameters characterizing the interaction between <p and x are small and/or in the 
decay of residual inflaton oscillations. In the particular case of parametric resonance 
modelled by a Mathieu equation, this contribution is indeed small. The growth of 
the modes of x is exponential in the resonance bands or tongues in parameter space, 
and the first resonance band is the only band wide enough to give rise to significant 
X excitations. However, this need not be so when there are more frequencies of 
excitation of the x field. Here we show that in this case resonance is governed by 
a general Hill equation and that other resonances beyond the first may contribute 
to the amplification of the x modes. Therefore, there is the possibility that these 
higher frequency excitations contribute significantly to the overall creation of x 
particles. This analysis is the subject of the present work, in the case when the x 
field is parametrically forced by two frequencies in a 1 : 2 ratio. For this case, we 
show that the contribution for the x field energy density of the higher harmonic is 
a small fraction of that of the fundamental frequency. 

The outline of the paper is as follows. In Section 2 we review the parametric fre- 
quency pre-heating mechanism and the method to compute the particle production 
rate in the general framework of Hill's equation. In Section 3 we apply this method 
to compare the reheating efficiency of two different couplings of the inflaton field 



to the x field. In Section 4 we sum up the conclusions of this analysis. 

2 Particle production by parametric resonance 

We start by reviewing some properties of the parametric resonance mechanism for a 
general periodic perturbation of the frequency of the oscillator, which corresponds 
to Hill's equation (J3j). 

Reheating models in inflationary universes start by considering that, at the end 
of inflation, the inflaton is in a coherent oscillatory state described by a space- 
independent expectation value, governed by the equation of motion for the inflaton 

4> + 3H<j> + a' 2 V 2 (j) + V ((/>) = 0, (1) 

where H is the Hubble parameter of a Friedman-Robertson- Walker metric with 
scale factor a(t), and V'(4>) is the derivative of the inflaton's potential V((f>) with 
respect to <f>. It is assumed that V{4>) has a vanishing minimum for the oscillations 
to take place, and also that the inflaton couples to another scalar field % which is 
then periodically perturbed by the inflaton. The equation of motion of x is given 
by 

X + SH X + a' 2 V 2 X + U'ix) + ^-V int (<i>, x) = 0, (2) 

OX 

where Vi n t = Vi n t(4>,x) is the interaction potential between \ and (f>, and where 
U(x) is a potential that gives mass to x- At the onset of the process of energy 
tranfer from 0, x is assumed to be at the vanishing minimum of this potential. The 
oscillations of (j> around the minimum of its potential are faster than the expansion 
rate of the universe, so it is meaningful in a relatively short time scale to work 
within the simplifying assumption that the actual spacetime can be approximated 
with a Minkowski metric. For reasonable choices of a single interaction potential 
Vi n t, namely, cubic interactions of the form g<ftx 2 i or quartic interactions h(p 2 x 2 , 
the resulting x equation ([2]) is that of an oscillator with a harmonically perturbed 
frequency. This yields a Mathieu type equation and provides a well-known mecha- 
nism for the parametric resonance of x an d for the exponential amplification of its 
particle number [T71 EZj . 

For more general coupling terms, in particular non-homogeneous couplings in- 
volving different powers of (ft, the perturbation equation for % is parametrically 
forced by a periodic function (for the quasi-periodic case see for instance [H]) and 
the equations of motion for the x modes are of the form (more details in Section [3|) 

Xk + (f4 + eF(t)) X k = 0, (3) 

where F(t) is a periodic function with period 2n, and we may set F(t) = F(—t), 
and L n F(t)dt = 0, and e a small positive parameter related with the amplitude of 
the inflaton oscillations. 

For the general Hill equation 0, Floquet's theorem states that the solutions 
are of the form 

Xk(t) = e^ X f\t), (4) 

where x k v) is a periodic function with the same period as F(t) and jik is one of 
the two characteristic exponents, which are both real or complex conjugate. 



Clearly, the lines Re(^fc) = divide the (wfe,e) parameter plane in unstable 
regions and stable ones, defining instability bands which become narrow close to 
the e = axis, producing what is often called a structure of tongues [40] . For a 
perturbation of period 2tt these tongues end on the e = axis at the points uj k that 
satisfy the parametric resonance condition uj k = n/2 (see Figure [T]-a). These modes 
are amplified by an arbitrarily small parametric forcing of period 2tt. 

To determine the characteristic exponent of the solutions of equation §3§ we 
follow Hill's method of solution. Given that the function F is periodic we write it 
as 

p(t) = \J2 c ^ ht > ( 5 ) 

k=— s 

where cq = (due to the parity of F) and s £ N|J{oo}. This (finite or infi- 
nite) Fourier series expansion of F(t), ©, together with Floquet's theorem suggest 
looking for a solution of ([3]) of the form 

Xk(t) = e^ Y, b n e M . (6) 

n=— oo 

Inserting this expression and equation ([5]) in ([3j), and equating the coefficients of 
e (fi k +tn)t^ we derive a homogeneous system with infinitely many linear equations 

(fi k + in) 2 b n + - ^ c m b n - m = 0, n = ... ,-2,-1,0, 1,2, ... (7) 

m=— oo 

Eliminating the coefficients b n in ([7]) , a non trivial solution exists if the characteristic 
exponent satisfies an infinite determinantal equation (called Hill's determinantal 
equation) , 

A(e,n k ) = \B rs \ = (8) 

where the elements B rs of A(e,/i^) are given by 

{1 if r = s 

(H + ir) 1 + u k 

Here an infinite determinant D = \B mn \, (m,n = — oo, . . . , +oo) is defined as the 
limit of D m = det(Bij)(i,j = —m, . . . , m) asra-> +oo, if it exists. 

Equation ([8]) can be reduced to the simpler form of a transcendental equation 
in/x fc @U 

sin 2 (i7r^fc) = A(e, 0) sin 2 (7rwfc). (10) 

This determines a characteristic exponent [i k , which in turn determines the b n 
coefficients of (J7|), and hence, a formal solution ([6]) of Hill's equation. The real part 
of /Xfc determines the growth factor of the solutions ([6]) . 

In the narrow resonance regime the phenomenon of parametrically resonant 
excitations, where e <C 1, it is possible to derive an approximate expression in closed 
form for \x k as a function of e. First, notice that A(0, 0) = 1. Second, we see from ([9]) 
that the values oiB rs depend linearly on e, and that A(e, 0) = l+0(e 2 ). This means 



A(e,0) ~ 1 + ae 2 , where a depends on c = (c_ s , c_ s +i, . . . ,c_i,ci, . . . , c s _i, c s ) and 
on a;^. 

Inverting (I10p using arcsin z = In ( iz + yl — z 2 J and expanding in powers of e, 
yields for the growth factor 

Re(fj, k )~ea(c,u} k )\sm(nu k )\. (11) 

This expression will be used later to determine numerically the resonant tongues 
for a particular function F(t). 

For periodic perturbations with finite Fourier expansions ©, s G N, it can be 
shown that the asymptotic form of the width L n of the re-th interval of instability 
is given by [H] 

L n = aJ-^- + 0{e p+1 ), s{p-l)<n<sp, (13) 

where p is the integer defined in these equations and where a n is a constant indepen- 
dent of e. One of the remarkable consequences of these formulae is that the widths 
of all the resonant bands up to the s-th band depend linearly on the perturbation 
amplitude, while the resonances of order higher than s are associated with thinner 
instability regions. 

This has direct impact in the computation of the energy production by para- 
metric resonance. For small e the main contributions to the energy density of the 
X field 

PxV) = \fd 3 k[u 2 k \ X k\ 2 + \Xk\ 2 ] (14) 

d^l\ Xk {t)\ 2 (15) 



come from the bands with widths depending linearly on the perturbation amplitude 
e, since the other have relatively negligible widths due their nonlinear cuspidal form. 
This yields 



Px (t) ~ 4tt ^ k 2 u, 2 k L m \ Xk (0)\ 2 ^ ht 

m=l 



(16) 

aj fe ~m/2 



where the value of Xfc(O) is evaluated at the center of the resonances bands. 

If all the modes start out with an amplitude |xfc(0)| then, by the virial theorem, 
or assuming that there is one x particle on each mode, 

|X*(0)| 2 = -, (17) 

and so 



p x {t) ~4vr ^k 2 u k Lie 



.pUkt 



(18) 



The growth of the x held modes persists in an expanding universe, if the time 
scale for the resonance is much shorter than the expanding time scale [14] . If we 
change the field to 

ipk = a-Xk, (19) 

and introduce conformal time 77, defined by dr\ = ma~ 1 dt, the equation (|2|) for ift^ 
becomes 

< + ( -^2 + eF( V )a 2 - 4-1 1>k = 0, (20) 

where the prime denotes the derivative with respect to 77. 

Assuming that m^> H, we can treat the expansion of the universe adiabatically, 
and thus, at any given time, its effect is a shift in the oscillatory frequency 



a" 



k z ^k z , (21) 

a 

and also an adiabatic increase of the amplitude of the driving force. In the ex- 
panding universe the parametric resonance analysis is only applicable if several 
conditions are satisfied: (i) e/m? <C 1 (perturbative regime) , (ii) the expansion 
of the universe can be neglected, that is H <C m, and (Hi) we also need the time 
scale on which the instable solution grows to be smaller than H, so that expansion 
is unimportant, H/m <C e, and finally, (iv) that the frequency does not redshift 
out of the resonance band in a time interval shorter than the amplification period 
m/(2/ife). This latter condition, for the z-band, can be written in the form 



m d 

— Wfc 

2fi k drj 



< L u (22) 



where //& and oJt are evaluated at the center of the resonance ii/i-band. 

The regime where these conditions can be easily satisfied is the period after 
inflation where naturally one has H/m <C 1 and where the equation of motion, in 
conformal time, for <f> = acf), is given by 

4" + 4(l-S) = 0, (23) 

with 

1 a" 



m 2 a 



< 1. (24) 



3 Instability pockets in higher order resonances 

We now apply the results of the previous section to the x dynamics in pre-heating, 
by considering two cases: (i) a cubic interaction term, g4>x 2 -> between the <f> and x 
fields, and (ii) the more general case of a cubic plus a quartic interaction, g<ftx 2 + 
hcjP'x 2 - The former case is one of the models that yield a Mathieu equation for 
the x modes, and is here only briefly considered for the purpose of illustrating the 
usual analysis and for comparison with our extended model (ii), in which the x field 
equation is of the form (j^J) with F(t) given by (jSJ) and 8 = 2, the simplest possible 
extension of the usual single frequency Mathieu model. According to (|13|) and (|18p 

6 



one would expect that, for comparable cubic and quartic terms, the efficiency of 
the energy transfer process should approximately double in the more general case. 
Assume that the background field (j) oscillates with frequency large compared 
with the Hubble expansion rate, which is always satisfied asymptotically \JA\ 143] . 
If we neglect the expansion in the eft dynamics then the solution of ([T]) can be given 
by, with a specific set of initial conditions, 



4>(t) = Acos(mt). 



(25) 



3.1 Single frequency interaction: the Mathieu equation 
case 

Consider the interaction potential g4>X 2 ■ This yields for the /c-mode equation, using 
r = mi, 

" k 2 gA 



xZ + 



„ H 9 cos(r) 

m m z 



Xk 



0. 



In this case one has s = 1, and, 



2 
w fc = 9 



e 



1 



m' 



(26) 

(27) 
(28) 
(29) 



and thus the width of the first instability band, when it is linear on e as given by 
(ED, is 



L 



_g± 

4m 2 



Using (I18p . the energy density in this case is given by 

7T , ,,0, 



aw = ;^, 



(30) 



(31) 



where fJ[ is the characteristic exponent of the first band, evaluated at the center of 
the band at height e. 



3.2 Multi-frequency interactions: the Hill equation case 

For the interaction potential geftx 2 + h(f) 2 x 2 the equation for the fe-th x mode is, 
with t = mt, 



xl + 



k 2 



hA 2 gA 



in 



2 2m 2 m 2 



hA 2 



H ^ cos(r) + - — k cos(2t 



2m 1 



In this case one has d3l) and d5l) with s = 2, and 



ul 



ci 

C'2 



m 2 2m 2 
1 

hA 

gA 



m 

7 



2' 



Xk 



0. 



(32) 

(33) 
(34) 
(35) 

(36) 



Equations (I13p for the widths of the first and second instability bands, are given 
by (|12p an (|13|) in the regions where they are linear in e, and yield 

L ' = & < 37 > 

i2 = is?- (38) 

The energy density in this case is given by, using (|18p . 

pi (r) = - 9 Ae^ r / m + -/M 2 e^ T / m , (39) 

where //5 and /^ are the characteristic exponents of the first and second bands, 
evaluated at the center of the bands at height e. 

3.3 Reheating efficiency of the two couplings 

To proceed further in the comparison of the growth of the \ field energy density 
generated by the two different couplings considered in this section, we must compute 
the relevant characteristic exponents in both cases. 

In order to determine the value of the real part of the characteristic value we 
use the equation (|lip . where the matrix [B rs ] was truncated at a size of 11 x 11 
after numerical accuracy tests. Then we have for the value of a as a function of Uk 
and of C2 



a(c 2 ,(jJk) 




n?){to 2 k -(n+lY) 

4 \ 



n?){ujl-{n + 2f 



(40) 



A(wl - l) 2 
and for the value for the growth factor 

Re(//fe) = easin(7ru;fc)- (41) 

Notice that if we set c 2 = in the latter expressions we recover the Mathieu case. 

In Figured] we show the stability diagrams for the two cases given above. Notice 
the appearance of an 'instability pocket' in the second resonance band of equation 
(j32[) in Figure [Hb). This phenomenon, which is much less well known than the 
appearance of additional non-cuspidal instability tongues associated with higher 
harmonics of the parametric forcing term, was studied in depth in [4HI42J. It plays 
a major role in explaining why the contributions of higher order resonances may be 
neglected and why the x field excitations are essentially single mode. 

As shown in Figure [H there are two major differences when comparing the 
bifurcation diagram of the two-frequencies Hill case with the Mathieu case. On the 
one hand, the width of the first instability band is slightly larger and its level curves 
of constant /i are slighted tilted when compared to the Mathieu case (Figure QJa). 
On the other hand, the shape of the second instability band is distorted giving rise 
to the emergence of a pocket, and its level curves of constant \i are lifted up. Also, 
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Figure 1: (a) Instability diagram for Mathieu equation <PS|) : (b) Instability diagram for Hill equation 
(j52f . Notice the emergence of an instability pocket in the second band in Figure [T] b) , and the modifica- 
tions in the level curves of the characteristic exponent \i. These plots were obtained using Hill's method 
of solution. 



the values of /i crossed by straight lines of fixed e in the first instability band, are 
significantly larger than those crossed in the second band, with the exception of the 
small region close to e = 0. 

In order to compare systematically the efficiency of the two couplings in trans- 
fering energy from the inflaton to the x field, we have computed the time that it 
takes to reach an e-fold increase of the total number of particles as a function of 
the coupling strength e, for C2 = 0.5. Instead of equations ([ST]) and (fHUj) . which are 
valid only in the limit of small e, we use the general equations (|18j) with Lj taken 
as the numerical value of the band widths for each model. 
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Figure 2: Time of one e-fold increase in the total number of particles of the \ held a s a function of the 
coupling strength e, for C2 = 0.5. 



In Figure 2 a), we see that the e-fold time is essentially determined by the 



contribution of the first instability band, except in the region of very small values 
of e. In this region, the multi-frequency coupling becomes more efficient than the 
Mathieu model because of the contribution of an additional linear instability band, 
but the effect has no cosmological implications since the energy transfer achieved 
by both couplings is negligible for these parameter values. 

In Figure 2.b) the overall behavior of the e-fold time as a function of e is shown 
for the two models. The kinks that can be seen in the two curves correspond to the 
values of e where the first band hits the u)k = axis. Due to the formation of the 
pocket in the second instability band, the multi-frequency model becomes actually 
less efficient than the single frequency excitation for moderate and large values of 
the coupling. 

This effect is also of limited cosmological relevance, since the e-fold times of the 
two models are of the same order of magnitude for similar values of the coupling 
strength. However, it is somehow unexpected that an additional linear (as opposed 
to cuspidal) instability band may translate into a less efficient resonance mechanism 
for most parameter values. This is of course a consequence of the pocket formation 
phenomenon, and it shows that the conclusions based on the analytic expressions 
for the asymptotic behavior of the instability bands for small values of e cannot be 
extrapolated. 

3.4 Reheating efficiency in the expanding universe 

We now take into consideration the expansion of the universe and discuss how this 
affects the contribution of the resonant bands found for the case of the two different 
couplings. 

The condition which should be satisfied for the resonance mechanism to work 
is given by equation (J22j) . As discussed in section [2] (see [H] ) , it translates the 
requirement that the frequency of the x particles should not be redshifted out of 
the resonance band in a time interval shorter than the amplification period m/(2^.). 

Using equation (f2Tj) . equation (f22j) can be recast as 



m 



l_(a" 

~2 



< K| ( ^) U (42) 

ml 



for the i-th instability band. 

According to the analysis of this section, significant particle production occurs, 
for either model, only for values of e such that the asymptotic approximations (II3p 
no longer hold, and that the 'instability pockets' of Figure l.b) are instead fully 
formed. Therefore, for a given e, the two-frequency model will exhibit a combination 
of lower values of //& and smaller width L2 than the single frequency model. Both 
effects contribute to making condition (I42|) harder to meet. For the same coupling 
strength and increasing the expansion rate, the drift across the second instability 
band will become swifter than the amplification time for the two-frequency model 
first. 

Therefore, when the expansion of the universe is taken into consideration, the 
conclusion that the reheating mechanism is weaker in the multi-frequency case than 
in the Mathieu case is reinforced. Since, however, only the contribution of the first 
instability band is accounted for in the efficiency estimates found in the literature, 
these estimates remain valid for the multi-frequency model. 
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Now equations (I3ip and (|39[) are valid as long as any given mode remains in the 
resonance band. Due to the expansion, the time interval At during which a mode 
remains in the band is 

At ^ < 43 » 

As long as the total time is small compared with H~ l /m, the total energy produced 
during the time interval NAt is approximately given by Np x . In this scenario, 
reheating is efficient if the ratio 

^ (44) 

P<p 

becomes of order one after a time smaller than the Hubble time. Otherwise a 
significant fraction of the original energy density is redshifted away. 

For the first band (as we have seen, the only one that might contribute), the 
latter quantity (|44j) reaches the value one for 



m 



A 2 



N = ^ ~ : . . (45) 

Px IgAe^l™ 

Therefore the condition for sufficient reheating NAt < H^ 1 /m becomes 

z}l e ^ T / m > J^ 2 > (46) 

or, 

^— > w (i^™ 2 ) ' ( 4? ) 

m \n- i a J 

where W is Lambert function [35] • Equations (|43|) and (|47|) establish a relation 
between the free parameters for sufficient reheating to take place that holds for 
both models considered in this section. 



4 Conclusions 

We have shown how the consideration of extra frequencies in the x equation of 
motion can be found in a simple model with V((j)) = m 2 (p 2 /2 for the inflaton 
potential and the interaction potential V(<f), x) = 50X 2 + h(j) 2 x 2 - In the narrow 
resonance regime the phenomenon of parametrically resonant excitations of the 
scalar field x by the inflaton's oscillations is then governed by a Hill equation, 
where the forcing term in the x equation has two distinct and commensurable 
frequencies. 

As a result of the presence of various harmonics in the equation of motion of the 
X field, the geometrical features of the resonant bands in the bifurcation parameter 
space are modified. Two main changes take place with respect to the single fre- 
quency (Mathieu) case. On the one hand, for small amplitudes of excitation, there 
are two resonant tongues with linear dependences on the amplitude, rather than 
just one as in the Mathieu case. On the other hand, closely related to the previous 
effect, there is a distortion of this additional band which gives rise to the formation 
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of 'instability pockets'. Due to the simultaneous presence (and interference) of two 
excitation frequencies, the lines in parameter space that correspond to the periodic 
solutions (/j, = 0) and define the boundary of the instability band cross each other 
for values of the forcing amplitude of order one. 

For the two cases under consideration we have evaluated and compared the 
particle production rates. We have considered first that the inflaton field 4> oscillates 
around the minimum of its potential much faster than the expansion rate of the 
universe, m^ 3> H, so that the expansion of the universe may be neglected. We 
have shown that, in general, the presence of an additional excitation frequency 
hinders, rather than favors, the efficiency of parametric resonance as an energy 
transfer mechanism, and that this is a consequence of the 'instability pockets' in 
the bifurcation diagram of the general Hill's equation. The enhancement of particle 
production due to the presence of a second linear, as opposed to cuspidal, instability 
band is shown to occur only for extremely small coupling strengths, for which 
both models yield negligible rates of \ particle creation. We then argue that the 
effects of an expanding universe further justify the approximation of neglecting the 
contribution of the second non-cuspidal instability band. 

In conclusion, our detailed analysis of two different coupling terms supports and 
justifies the usual approach in the literature, where the efficiency of reheating by 
parametric resonance is evaluated by considering the simplest form of the paramet- 
rically forced equation, and only the dominant contribution of the first instability 
band. 
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